Wave modelling in a cylindrical non-uniform helicon discharge 
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A radio frequency (RF) field solver based on Maxwell's equations and a cold plasma dielectric tensor is em- 
ployed to describe wave phenomena observed in a cylindrical non-uniform helicon discharge. The experiment 
is carried out on a recently built linear plasma-material interaction machine: the MAGnetized Plasma In- 
teraction Experiment (MAGPIE) [B. D. Blackwell, J. F. Caneses, C. Samuell, J. Wach, J. Howard, and C. 
S. Corr, submitted on 25 March 2012 to Plasma Sources Science and Technology], in which both plasma 
density and static magnetic field are functions of axial position. The field strength increases by a factor of 
15 from source to target plate, and plasma density and electron temperature are radially non-uniform. With 
an enhancement factor of 9.5 to the electron-ion Coulomb collision frequency, 12% reduction in the antenna 
radius, and the same other conditions as employed in the experiment, the solver produces axial and radial 
profiles of wave amplitude and phase that are consistent with measurements. Ion-acoustic turbulence, which 
can happen if electron drift velocity exceeds the speed of sound in magnetized plasmas, may account for the 
factor of 9.5 used to match simulated results with experimental data. To overcome the single m vacuum solu- 
tion limitations of the RF solver, which can only compute the glass response to the same mode number of the 
antenna, we have adjusted the antenna radius to match the wave field strength in the plasma. A numerical 
study on the effects of axial gradient in plasma density and static magnetic field on wave propagations is 
performed, revealing that the helicon wave has weaker attenuation away from the antenna in a focused field 
compared to a uniform field. This may be consistent with observations of increased ionization efficiency and 
plasma production in a non-uniform field. We find that the relationship between plasma density, static mag- 
netic field strength and axial wavelength agrees well with a simple theory developed previously. A numerical 
scan of the enhancement factor to the electron-ion Coulomb collision frequency from 1 to 15 shows that the 
wave amplitude is lowered and the power deposited into the core plasma decreases as the enhancement factor 
increases, possibly due to the stronger edge heating for higher collision frequencies. 



I. INTRODUCTION 

Generically, a helicon discharge usually refers to a 
cylindrical plasma discharge with an axial static mag- 
netic field, driven by radio frequency (RF) waves at fre- 
quencies between the ion and electron cyclotron frequen- 
cies, ui C i <C u <C c^ce- 1 A helicon discharge produces plas- 
mas with densities typically much higher than capacitive 
and inductive plasma sources operating at similar pres- 
sures and input RF powers^ Because of this high ioniza- 
tion efficiency, helicon discharges have found applications 
in various fields, including: plasma rocket propulsion^— 
a plasma source for magnetic fusion studies^ Alfven wave 
propagation^ RF current drive^ laser plasma sources^ 
possibly semiconductor processing, electrodeless beam 
sources, and laser accelerators^ 

To date, most helicon studies have treated devices 
with uniform static magnetic fields, however, many ap- 
plications require operation with axial magnetic field 
variations.— A few researchers have investigated helicon 
plasma sources with non-uniform magnetic fields, and 
have found that the plasma density increased when a 
cusp or non-uniform magnetic field was placed in the 
vicinity of the helicon antennaJ^r— However, detailed 
examination of the reasons for this enhanced plasma den- 
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sity has not yet been conducted, although fast electrons 
and improved confinement are mentioned as possible con- 
tributors. Guo et aljii furthered this study by looking 
at the effects of non-uniform magnetic field on source 
operations, and found that strong axial gradient in den- 
sity associated with non-uniform field configuration can 
contribute to the absorption of wave fields and a high 
ionization efficiency. Takechi et alJ£ also suggested that 
there may be a close relationship between plasma density 
profile and RF wave propagation and absorption regions, 
finding the density uniformity in the radial direction im- 
proved markedly with the cusp field. Therefore, studying 
the effects of various static magnetic field configurations 
on helicon wave propagation is of significant importance 
to producing desired plasma profiles and understanding 
the role of magnetic field in helicon plasma generations. 

This paper is dedicated to modelling the wave field 
observed in MAGPIE (MAGnetized Plasma Interaction 
Experiment), and investigating helicon wave propagation 
in the non-uniform magnetized plasma of this machine, 
in which both the static magnetic field and its associ- 
ated plasma density are functions of axial position. The 
plasma density and electron temperature are also depen- 
dent on radius. We assume in this study that the electron 
temperature is independent of z and the static magnetic 
field is almost independent of r (Eq. (7)). MAGPIE is 
a linear plasma-material interaction machine which was 
recently built in the Plasma Research Laboratory at the 
Australian National University, and designed for study- 
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ing basic plasma phenomena, testing materials in near- 
fusion conditions, and developing potential diagnostics 
applicable for the edge regions of a fusion reactor.— A 
RF field solver, — based on Maxwell's equations and a 
cold plasma dielectric tensor, is employed in this study. 
The motivations of our work are to explain the wave 
field measurements in MAGPIE, and to study the effects 
of magnetic field configuration on helicon wave propa- 
gation. The rest of the paper is organized as follows: 
Section II describes the experimental apparatus and di- 
agnostic tools, together with the measured static mag- 
netic field, plasma density and temperature profiles; Sec- 
tion III provides an overview of the employed theoreti- 
cal model and the numerical code, together with com- 
parisons between computed and measured wave fields; 
Section IV is dedicated to a numerical study of the ef- 
fects of plasma density and static magnetic field profiles 
on the wave propagation characteristics; Section V aims 
to study the physics meaning of the enhancement fac- 
tor to electron-ion Coulomb collision frequency, and the 
effects of the direction of static magnetic field on wave 
propagations^ - — finally, Section VI presents concluding 
remarks and future work for continuing research. 



and 0.08 m in radius. Gases are fed through the down- 
stream end of the target chamber, and drawn to the up- 
stream end of the source tube by a 170 L/s turbo pump. 
Gas pressures are measured in the target chamber by a 
hot cathode Bayard- Alpert Ionization gauge (< 0.01 Pa), 
a Baratron pressure gauge (0.01-10 Pa) and a Convec- 
tron (0.1 Pa-101.33 kPa) for pressure process. In this 
experiment, argon gas is used with a filling pressure of 
Pb = 0.41 Pa. The two regions, source and target, are 
surrounded by a set of water cooled solenoids, with in- 
ternal radius of 0.15 m. These source and target sets of 
solenoids are powered by two independent 1000 A, 20 V 
DC power supplies, providing flexibility in the axial con- 
figuration of the static magnetic field, e. g. maximum 
of 0.09 T and 0.19 T in the source and target regions, 
respectively. The non-uniform field configuration is ex- 
pected to provide a flexible degree of radial confinement, 
better plasma transport from the source tube to the tar- 
get chamber, and possible increased plasma density ac- 
cording to previous studies* ^ 17 ' 19 The direction of the 
static magnetic field employed in the present work points 
from target to source. 



II. EXPERIMENT 

A. Experimental setup 

Similar to other helicon devices, 1 MAGPIE mainly 
consists of a dielectric glass tube surrounded by an an- 
tenna, a vacuum pumping system, and a gas feeding sys- 
tem, together with a power supply system connected to 
the antenna, and various diagnostics. Figure 1 shows a 
schematic and introduces a cylindrical (r, 9, z) coordi- 
nate system^ The plasma is formed in the region under 
the antenna (—0.243 < z < —0.03 m) and the near field 
to the antenna^ Following convention, however, we de- 
fine the whole glass tube (— 1 < z < m) as the source 
region and the compressed field region (0 < z < 0.7 m) 
as the target region (or equivalently "diffusion region" in 
some references). In MAGPIE, the z < —0.243 m region 
is named "upstream" and z > —0.03 m "downstream". 

A glass tube of length 1 m and radius 0.05 m is used 
to contain source plasmas in MAGPIE. A left hand half- 
turn helical antenna, 0.213 m in length and 0.06 m in 
radius, is wrapped around the tube and connected to a 
tuning box which can be adjusted between 7 and 28 MHz, 
a directional coupler, a 5 kW RF amplifier, and a 150 W 
pre-amplifying unit. For the present study, an RF power 
of 2.1 kW, 13.56 MHz, pulse width of 1.5 ms and duty cir- 
cle of 1.5 % is used. The antenna current is measured by 
a Rogowski-coil-type current monitor. For these experi- 
ments an antenna current of magnitude I a — 38.8 A was 
measured. A grounded stainless steel cylindrical mesh 
surrounding the whole source region is employed to pro- 
tect users. The source region is connected on-axis to 
the aluminium target chamber which is 0.7 m in length 



B. Plasma profile diagnostics 

A passively compensated Langmuir probe was em- 
ployed in our experiment to measure the plasma density 
and electron temperature, calculated from the I{V) curve 
obtained by an Impedans Data Acquisition system. 24 
The probe comprises a platinum wire of diameter 0.1 mm, 
and a surrounding alumina insulator. The length of the 
insulator is 6 mm shorter than that of the platinum wire 
so that the exposed platinum wire forms the probe tip. 
Electron currents were drawn to clean the probe during 
regular intervals of argon discharges. The probe is lo- 
cated at z — 0.17 m as shown in Fig. 1. 

Typical measured axial profile of field strength and ra- 
dial profiles of plasma density and electron temperature 
in MAGPIE are shown in Fig. 2. Particularly, as shown 
in Fig. 2(a), the increase in field strength (Bq(z)) from 
antenna end (z = —0.243 m) to field peak (z — 0.51 m) is 
a factor of 15. The axial profile of plasma density (n e (z)) 
is assumed to be proportional to Bq(z), consistent with 
generally accepted knowledge that the density follows the 
magnetic field linearly 25 i 26 Figure 2(b) shows the radial 
profiles of plasma density (n e (r)) and electron temper- 
ature (T e (r)), measured at z = 0.17 m and fitted with 
straight lines. During the density fitting procedure, in or- 
der to avoid negative fitted values, the density was set to 
zero in the region of 0.066 < r < 0.08 m. We assume the 
total density profile is separable, n e (r, z) = n e {r) x n e (z). 
The fitted lines in n e (r) and T e (r), and the measured 
Bq (z) data will be used in section III to constrain wave 
field simulations. 
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FIG. 1. A schematic of the MAGPIE (MAGnetized Plasma Interection Experiment).— A circle denotes the position of the 
helicon antenna which is left hand half-turn helical. The dot-dashed line is the machine and coordinate system axis, defining 
r = m. The coordinate system is right-handed with 8 — chosen to be the zenith angle. 
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FIG. 2. Typical measured profiles: (a) axial profile of static 
magnetic field on axis, (b) radial profiles of plasma density 
(dots) and electron temperature (squares) at z = 0.17 m, 
together with their fitted lines, solid and dashed, respectively. 
The solid bar in (a) denotes the antenna location. 



C. Wave field diagnostics 



Helicon wave fields were measured by a 2-axis "B dot" 
or Mirnov probe. Details about the probe can be found 
in Blackwell et al..— To measure the axial profiles of B r 
and B Zl the probe was inserted on axis from the end of 
the target chamber. The probe is long enough to mea- 
sure B r and B z in the range —0.25 < z < 0.7 m. Two 
perpendicular magnetic field components (B r and B z in 
this case) can be sampled simultaneously. To measure 
the radial profiles of the three magnetic wave compo- 
nents, B ri Bg and B z , the probe was inserted radially 
at z = 0.17 m, and rotated about its axis to measure 
Bg and B z . The B-dot probe couples inductively to the 
magnetic components of the helicon wave and electrostat- 
ically to the RF time varying plasma potential. To limit 
our measurements to the inductively coupled response, 
a current balun was employed to screen the electrostatic 
response. Further information about the procedure to 
eliminate the electrostatic response of the probe can be 
found in Franck et al.. 2 ? Both axial and radial profiles of 
wave phase were measured through a phase-comparison 
method, similar to Light et al.i2& To measure the vari- 
ation in wave phase with axial position, the signal from 
an on- axis axially inserted probe was compared to the 
phase of the antenna current. A similar procedure was 
conducted to measure the variation in wave phase with 
radial position at z = 0.17 m. It should be noted that all 
probe diagnostics are intrusive, and can affect the plasma 
parameters, and hence the wave fields. 
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III. SIMULATION 

An RF field solver (or ElectroMagnetic Solver, EMS) 20 
based on Maxwell's equations and a cold plasma dielec- 
tric tensor is employed in this study to interpret the RF 
waves measured in MAGPIE. This solver has been used 
successfully in explaining wave phenomena in two other 
machines: a helicon discharge machine at The Univer- 
sity of Texas at Austin^ 2 - and the LArge Plasma Device 
(LAPD) at the University of California at Los Angles^ 
Details of the solver can be found in Chen et al.,— while 
a brief overview is given below. 



A. Theoretical model 

The Maxwell's equations that this solver employs to 
determine the RF wave field in a helicon discharge are 
written in the frequency domain 



V xE = iujB, 



Mo 



V xB 



-iujD 



(1) 



(2) 



where E and B are the electric and magnetic fields, re- 
spectively, D is the electric displacement vector, uj is the 
antenna driving frequency, and j a is the antenna current 
density. The quantities D and E are linked to each other 
by a dielectric tensor e a p that represents vacuum, glass 
and plasma. In the vacuum and glass regions, the dielec- 
tric tensor is e a p = e*(r, z)6 a p, where 8 a p is the Kro- 
necker symbol and £*(r, z) is a scalar. The term e*(r, z) 
equals to 1 and e g for vacuum and glass regions, respec- 
tively, where e g is the dielectric constant of glass. In the 
plasma region, because of the cold plasma approxima- 
tion made here, the relation between D and E is in form 

o f29,30 



D = e (eE + ig\E x b] + (rj - e)(E ■ b)b), (3) 

where b = Bo/£>o is the unit vector along the static 
magnetic field and 



£ = 1 



v— » UJ + iv a 



uj (u) + iv a ) 2 — cj 



2 _ ,.,2 ' 
ca 



(4) 



Here the subscript a la bels particle species, i. e. elec- 
tron and ion, uj pa = \J n a q\l Enm a is the plasma fre- 
quency, uj ca = q a Bo/m a gyrofrequency, and v a collision 
frequency between species. The plasma is assumed to 
be nearly fully ionized in the present study, so that neu- 
tral collisions are neglected. Because u ee and va do not 
contribute to the momentum exchange between electron 
and ion fluids^ collision frequencies for electrons and 
ions species are v e = v e i — 2.91 x 10~ 12 n e T e 3 ^ 2 lnA and 
v i = v ie. = rn e m ^ ll/ ei, respectively, from which we can 
see Vi e v e i- Here, T e and n e are given in eV and m~ 3 , 
respectively, and the Coulomb logarithm is calculated to 
be In A = 12. Singly ionized argon ions are assumed in 
this study, so that qi = —q e = \e\. 

The externally applied Bo(r,9,z) is assumed to be 
axisymmetric, with Bo r <C Bq z and Bog — 0. There- 
fore, it is appropriate to use a near axis expansion 2 ^ for 
Bq(t, 9, z), namely Bq z is only dependent on z and 



B 0r (r, z) 



1 dB 0z (z) 

--T . 

2 dz 



(7) 



The antenna, as described in section lll Al is a left hand 
half-turn helical antenna. We assume that the antenna 
current is divergence free, to eliminate the capacitive cou- 
pling. Fourier components of the antenna current density 
are given by 



jar 0; 



(8) 



ja9 = la 1 



~ 5 ( r - R a)i±;{Kz - Za) + 5(Z - Z a - L a ) 



H(z-z a )H(z a +L a ~z) -imn[l+(z-Za)/Lg] 



), 



i — I - 

Jaz — ± a 



-5(r-R a ) 



(9) 



(10) 



xH(z - z a )H(z a + L a 



z 



where L a is the antenna length, R a the antenna radius, 
z a the distance between the antenna and the endplate 
in the source region, and H the Heaviside step function. 
Note that the antenna geometry selects only odd har- 
monic mode number to, as indicated by Chen ct al.j^ - 



B. Boundary conditions 



rj = 1 



pa 



uj {uj + iv a ) 2 - uj 2 ca ' 



uj(uj + iv a ) 



(5) For a given azimuthal mode number to, Eq. (1) and 

Eq. (2) are first Fourier transformed with respect to 
the azimuthal angle, and then solved through a finite 
difference scheme on a 2D domain (z, r), as shown in 

(g) Fig. 3. In the experiment, there is a radial air gap ( 
0.055 < r < 0.0585 m) between the antenna and the 
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glass tube, which is taken as glass region in the compu- 
tational domain. We found that simulated results are 
insensitive to the dielectric constant in the glass region 
0.05 < r < 0.055 m by varying the constant from 1 to 10 
and no change detected in the wave field, therefore, we 
expended the glass area radially to fill this air gap. The 
thickness of the antenna is approximately 0.002—0.003 m. 
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FIG. 3. Computational domain employed to simulate the ex- 
perimental setup shown in Fig. 1. Here all dimensions are 
given in millimetres. 
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The radial wall of the target chamber and the axial 
endplates are ideally conducting so that the tangential 
components of E vanish at the surface of these bound- 
aries, i.e., 



Eg(L r , z) = E z (L r ,z) = 0, 



(11) 



FIG. 4. Variations of magnetic wave field in axial direction 
(on-axis): (a) |_Br|rms, (b) phase of B r . Computed results 
(lines: dotted for v e g — v e i and Raim = Rc* P , dashed for 
f c ff = and -Rsim = 0.88i?cx P , and solid for z/ c ff = 9.5v e i 



and 7?aim 
(dots). 



0.88-R C x P ) are compared with experimental data 



E r (r,0) =E {r,O) = 0, (12) 



E r (r,L z )=E e (r,L z ) = 0, (13) 

where L r and L z are the radius of the target chamber 
and the length of the whole machine, respectively. More- 
over, all field components must be regular on axis, thus, 
Bg\ r= o — and (rEg)\ r= o — for m — 0; E z \ r= o and 
(rEg)\ r —o for m ^ 0. 21 In the present work, we choose 
the fundamental odd mode number m = 1, which is pref- 
erentially excited in the helicon discharge launched by a 
left hand half-turn helical antennai 23 i 28 ' 32 



C. Computed and measured wave fields 

Based on the measured field strength configuration and 
plasma profiles shown in Fig. 2, simulations are per- 
formed. Figure 4 shows the axial profiles of the computed 
B r amplitude and phase on axis, and their comparisons 
with experimental data. With the collisionality set to 
Z4ff = v eii where v e ft is the effective collision frequency, 
and the antenna radius set to match the experiment, the 
predicted wave field is ~ 30% of the measured value, 
and the profile a poor match to the experiment. It is 



possible to obtain better agreement by varying the colli- 
sionality, which strongly affects the profile but leaves the 
magnitude largely unchanged, and the antenna radius, 
which strongly affects the field amplitude and leaves the 
radial and axial profiles of B unchanged. A qualitative 
match between measurement and simulation of the axial 
variation of B r is found using an enhancement in colli- 
sionality of f e ff = Ci^ei + v- le ) ss (V e i with £ = 9.5, and 
an adjustment in antenna dimension of i? S i m = £i? cxp 
with £ = 0.88. Calculation of the axial gradient of the 
computed phase variation shows a travelling wave, with 
a good agreement with data. 

A number of physics reasons, detailed in Section [V] 
exist to support an enhancement in the collisionality. 
However, experimental uncertainties in the current, an- 
tenna and vessel dimensions are not sufficient to explain 
the 12% reduction in antenna radius required to find a 
qualitative match in field amplitude. A possible omis- 
sion of the RF solver, that might explain the need for 
an artificially reduced antenna dimension is the single m 
vacuum solution limitation. Specifically, EMS can only 
compute the glass response to the same mode number of 
the antenna. In reality, however, the antenna will gen- 
erate to ^ 1 vacuum harmonics. It is possible that the 
evanescent length scales of m = —1 and m = are much 
larger than 777 = 1, producing substantial 777 = — 1 and 
777 = waves at the plasma-glass surface. Mode coupling 
at the plasma-glass surface may then couple 777 =/= 1 to 
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m = 1, producing an increased field in the plasma. 

The local minimum observed both experimentally 
and numerically in the axial profile of |-B r | rms around 
2 = 0.27 m has been also observed in many other 
devices , 12 i 17 ' 28 i 33 for both uniform and non- uniform field 
cases. For the uniform field case, it has been sug- 
gested that the spatial modulation of the helicon wave 
amplitude is not caused by reflections from the end 
boundaries, but by a simultaneous excitation of two ra- 
dial modesi n ' 12 i 23 ' 28 Similarly, the minimum observed in 
MAGPIE cannot be explained by standing waves, be- 
cause the amplitude becomes much smaller at bigger z 
(suggesting strong damping), and the phase advances 
with increasing z (denoting a travelling wave). Further, 
radial profiles of the wave field in Fig. 5 feature a possible 
superposition of the first and second radial modes of the 
to = 1 azimuthal mode. Therefore, we speculate that the 
minimum observed here may be also due to the simulta- 
neous excitation of two fundamental radial modes. We 
will show later that radial gradient in plasma density is 
essential for the excitation of this local minimum under 
the present experimental conditions. 

Figure 5 shows the radial profiles of computed wave 
fields for v c r — 9.5v e i and i? S i m = 0.88i? CX p at three axial 
positions in the target region, together with the experi- 
mental data measured at z = 0.17 m. The predicted wave 
field amplitude profile at z = 0.17 m is consistent with 
the data, but the magnitude is nearly double the mea- 
sured value, and the phase profile is a poor match. We 
have also computed the wave fields at axial locations with 
best agreement to the amplitude (z = 0.21 m) and phase 
(z = 0.07 m). We justify this freedom of choice by the 
experimental uncertainty in axial density profile, and the 
numerical sensitivity identified in the radial profile of the 
wave field with axial position. Inspection of Fig. 5 reveals 
that it is possible to find a reasonable agreement to the 
wave amplitude and phase profile-albeit-independently. 
As expected, all calculations show the wave mode struc- 
ture of to = 1 through \B z (r = 0)| rms « 0, consistent 
with the antenna parity. The near-null minima in both 
\Bg(r)\ rms and |B z (r)| rms suggest a likely simultaneous 
excitation of the first and second radial modes. This 
may account for the minimum in Fig. 4(a), and is con- 
sistent with conclusions from othersj 1:L ' 12 ' 23 ' 28 i 32 Mori et 
alii 2 - suggested that the superposition feature is associ- 
ated with wave focusing effect caused by the non-uniform 
magnetic field in the target region. 



IV. NUMERICAL PROFILE SCANS 

It has previously been shown that the plasma density 
can be further increased by introducing a cusp or non- 
uniform static magnetic field in the vicinity of the helicon 
antenna?^— To shed light on the increased plasma pro- 
duction, we perform a detailed numerical study on the 
effects of radial and axial plasma density gradients and 
axial magnetic field gradient on wave propagation char- 



acteristics. The enhancement of £ = 9.5 to v e i and the 
adjustment of £ = 0.88 to R exp are still employed in this 
section because they provide a good agreement with the 
measured wave field. 



A. Axial profile of plasma density 

We first study the effect of axial gradient in plasma 
density, which has been assumed to be linear with the 
static magnetic field so far, on wave propagations by com- 
paring the wave fields from three different on-axis density 
profiles shown in Fig. 6. Other conditions are kept the 
same as previous sections. 

The computed wave fields in axial direction (on-axis) 
are shown in Fig. 7. A log scale in the amplitude has been 
employed to see the wave propagation features clearly. 
We can see from the phase variations (Fig. 7(b)) that as 
density is decreased in the target region the wavelength 
increases, which is consistent with a simple theory devel- 
oped previously, 11 

3-83 _ to n e e^ Q 
Rp k Bg 

Thus, if u), R p (plasma radius) and Bq are all fixed, k 
behaves proportional to n e , which means that the wave- 
length becomes larger at lower density. Here, the value 
of 3.83 is the first non-zero Bessel root of J\{r) — 0, rep- 
resenting the first radial mode, which is assumed to be 
dominant in our case. Further, for all density profiles 
shown in Fig. 6, the wavelength is bigger in region of 
< z < 0.6 m than that in other regions, indicating an 
increased phase velocity. This increased phase velocity 
together with the strong decay of wave amplitude sug- 
gest strong coupling of RF power from the antenna into 
the plasma in this regionJi 



B. Axial profile of static magnetic field 

Second, following section lTV A[ we keep the axially uni- 
form density profile and study the effects of axial gradi- 
ent in static magnetic field, which is radially near uni- 
form according to Eq. (7). Three employed field profiles 
are shown in Fig. 8, which enable us to study the ef- 
fects of field gradient in target and source regions seper- 
ately on wave propagations. Comparison between solid 
and dashed lines in Fig. 9 shows that axial gradient in 
magnetic field in the target region increases the prop- 
agation distance of helicon waves, consistent with Mori 
et al.'s conclusion that a focused non- uniform magnetic 
field provides easier access for helicon wave propagations 
than a uniform field^ The simple theory in Eq. (14) is 
satisfied again here: with decreased field strength in the 
target region, the wavelength becomes shorter. Although 
the difference between dashed and dotted field profiles is 
small as shown in Fig. 8, the computed field amplitudes 
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FIG. 5. Variations of magnetic wave field in radial direction: (a), (c) and (e) are |S r | rms , |i?e| rms and |-B z | rms , respectively; 
(b), (d) and (f) are the corresponding phase variations. Dots are experimental data while lines (dotted: z = 0.07 m, dashed: 
2 = 0.17 m, solid: z = 0.21 m) are simulated results. 
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FIG. 6. Normalized axial profiles of plasma density on axis. 
The solid line is associated with density profile linear with 
Bo(z). The dashed line is the same to the solid one except 
in region of < z < 0.7 m where the density maximum 
is adjusted to match preliminary experimental observations. 
The dotted line represents a z-independent density profile. 
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are significantly different. Figure 9 shows that with the 
uniform field profile wave amplitude is much bigger than 
that with non-uniform field profile for z < m. Fur- 
thermore, waves keep their travelling features till the left 
cndplatc for uniform Bq, whereas for non- uniform Bq(z) 
the wavelength becomes smaller when approaching left, 
and the waves are not travelling at all when Bo(z) is low 
enough (z < —0.9 m). 



FIG. 7. Variations of magnetic wave field in axial direction 
(on-axis) with n e axial profile: (a), log scale of |-B r | r m S ; (b), 
phase of B r . The three density profiles employed correspond 
to those in Fig. 6. 



C. Radial profile of plasma density 

Now, we keep the plasma density and static magnetic 
field both uniform in the axial direction, and study the 
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FIG. 8. Axial profiles of static magnetic field. The solid 
line shows original experimental data (Fig. 2(a)), while the 
dashed line shows the same except being flattened in region 
of z > m and the dotted line is flattened everywhere. 




FIG. 9. Variations of magnetic wave field in axial direction 
(on-axis) with Bo axial profile: (a), log scale of |-B r | r m S ; (b), 
phase of B r . The three field profiles employed correspond to 
those in Fig. 8. 




FIG. 10. Radial profiles of plasma density. 




z(m) 



FIG. 11. Variations of magnetic wave field in axial direction 
(on-axis) with n e radial profile: (a), log scale of |_B r |rms; (b), 
phase of B r . The two density profiles employed correspond 
to those in Fig. 10. 



V. COLLISIONALITY AND FIELD DIRECTION 



effects of radial gradient in plasma density. The two den- 
sity profiles employed are shown in Fig. 10, with and 
without radial gradient, and their corresponding results 
are shown in Fig. 11. We can see first that the local min- 
imum in wave amplitude profiles, e. g. at z = —0.52 m 
and z = 0.27 m, disappear when the radial density pro- 
file is fiat, suggesting that the radial gradient in plasma 
density is essential to have a local minimum under the 
present conditions. Second, the wave amplitude is much 
bigger in both target and source regions for plasma den- 
sity with radial gradient, suggesting that a radial gra- 
dient in density may be useful to maximize the plasma 
production. 



A. Enhancement of electron-ion collision frequency 

In a similar manner to other work) 21 i 22 we have used an 
enhancement to v e i (here v c r rj 9.5^ e i), in order to find 
a qualitative match of simulated wave field to the data. 
In this section, we explore the physical consequences of 
scaling f e ff/ v ei in simulations while keep using the ad- 
justment £ = 0.88 in the antenna radius. 

Variations of wave amplitude on axis in the axial di- 
rection for different collision frequencies are shown in 
Fig. 12. As the collision frequency increased from v e i to 
15z/ e i, the wave amplitude decreases nearly everywhere, 
and the wave decay length is shortened. This indicates 
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FIG. 12. Variations of on-axis wave amplitudes in axial direc- 
tion with different electron-ion collision frequencies: log scale 
of |B r | 

rms ■ 

that the wave energy or power coupled from the antenna 
to the core plasma drops as the collision frequency be- 
comes higher, and the power is more absorbed under the 
antenna. This is consistent with a previous conclusion 
that the RF energy is almost all absorbed in the near re- 
gion of the antenna rather than in the far region^ The 
oscillations near the downstream and upstream ends at 
low v e f{ are caused by reflections from the ideally con- 
ducting endplates, which disappear if the endplates are 
moved further away. 

As suggested by Lee et al.^ an enhanced electron- 
ion collision frequency may be due to ion-acoustic tur- 
bulence which can happen if the electron drift veloc- 
ity exceeds the speed of sound in magnetized plasmas. 
Based on the experimental conditions in MAGPIE, we 
have calculated the threshold field strength Bt, below 
which ion-acoustic turbulence can happen. This thresh- 
old is given by vd > C s , where vd « kBT e /\e\BoRp is 
the electron drift velocity with ks Boltzmann's constant 
and C s = ^ksTe/rrii the speed of sound in magnetized 
plasmas, resulting in Bt < 0.0224 T. Thus, the whole 
source region which produces helicon plasmas and waves 
is indeed located within this range. The ion-acoustic tur- 
bulence has the effect of providing additional electron-ion 
collisions within a dielectric tensor model, and thereby 
improves the agreement with observations. 



B. Direction of static magnetic field 

Observations have been made previously that the di- 
rectionality of helicon wave propagations is dependent 
on the direction of static magnetic field in helicon dis- 
charges using helical antennas, but all in uniform field 
configurations ! n i 22 i 23 i 34 l n this section, we study the di- 
rectionality in a non-uniform field configuration. Specif- 
ically, we have computed the wave amplitude and wave 
energy density in MAGPIE for the experimental and field 
reversed configurations. In MAGPIE, the field points 
from target to source, as mentioned in section lll Al Figure 
13 shows the computed axial profiles of wave amplitude 
on axis and 2D contour plots of wave energy density for 



both field direction pointing from target to source (Fig. 
13(a) and 13(b)) and field direction pointing from source 
to target (Fig. 13(c) and 13(d)). In this calculation, we 
have chosen v c g = v e n to see more details, and chosen 
the density profile to be linear with Bq(z) in the axial 
direction and non- uniform in radius as measured in Fig. 
2(b). The field strength profile used here is shown in Fig. 
2(a). Inspection of Fig. 13 reveals that the wave energy 
is larger in the opposite side of the antenna, relative to 
the direction of the static magnetic field. This observa- 
tion has been confirmed experimentally through finding 
that the plasma is brighter in the opposite side of the 
antenna relative to the direction of the applied external 
field. Therefore, the dependence of the direction of heli- 
con wave propagations to that of static magnetic field still 
exists even when the field configuration is non-uniform. 



VI. CONCLUSIONS 

A RF field solver based on Maxwell's equations and a 
cold plasma dielectric tensor is employed to describe the 
wave phenomena observed in a cylindrical non-uniform 
helicon discharge, MAGPIE. Here, the non-uniformity is 
both radial and axial: the plasma density is dependent 
on r and z, the static magnetic field varies with z, and 
the electron temperature is a function of r. A linear fit- 
ting was conducted for radial profiles of plasma density 
and electron temperature, and the fitted profiles utilized 
in wave field calculations. A linear relationship between 
the axial profile of plasma density and the static magnetic 
field was assumed. Other conditions used in the simula- 
tion were taken from experiment directly, including filling 
gas (argon), antenna current 38.8 A (magnitude), driv- 
ing frequency 13.56 MHz, and a left hand half-turn helical 
antenna. 

With an enhancement factor of 9.5 to the electron- 
ion Coulomb collision frequency v e i and 12% reduction 
in the antenna radius, the wave solver produced consis- 
tent wave fields compared to experimental data, includ- 
ing the axial and radial profiles of wave amplitude and 
phase. Particularly, a local minimum in the axial pro- 
files of wave amplitude was observed both experimentally 
and numerically, agreeing with previous studieSi 12 ' 17 ' 28 i 33 
Mode structure of m = 1 is consistent with the left hand 
half-turn helical antenna being used. A possible expla- 
nation for the enhanced electron-ion collision frequency 
has been offered through ion-acoustic turbulence, which 
can happen if electron drift velocity exceeds the speed 
of sound in magnetized plasmas. 22 Through calculating 
these two speeds based on MAGPIE conditions, we found 
that it is indeed satisfied in the source region of MAG- 
PIE where the helicon plasmas and waves are produced. 
Furthermore, to overcome the single m vacuum solution 
limitations of the RF solver, which can only compute the 
glass response to the same mode number of the antenna, 
we have adjusted the antenna radius to match the wave 
field strength in the plasma. 
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FIG. 13. Axial profiles of magnetic wave field (on-axis) and contour plot of wave energy density in (z, r) space for a non-uniform 
plasma density: (a), log scale of |-Br|rma for upstream Bq{z); (b), wave energy density for upstream Bo(z); (c), log scale of 
|£>r|rms for downstream Bo(z); (d) wave energy density for downstream Bq(z). 



A numerical study on the effects of axial gradients in 
plasma density and static magnetic field on wave prop- 
agations was carried out. This showed that the axial 
gradient in magnetic field increases the decay length of 
helicon waves in the target region. The strong decay of 
wave amplitude and the increase in phase velocity in re- 
gion of < z < 0.6 m indicate strong coupling of RF 
power from the antenna to the plasma, which is consis- 
tent with a previous study by Guo et al.J^ The relation- 
ship between plasma density, static magnetic field and 
axial wavelength is consistent with a simple theory de- 
veloped previously,!! 

A numerical scan of the enhancement factor to v e i re- 
veals that with increased electron-ion collision frequency 
the wave amplitude is lowered and more focused near the 
antenna. This is mainly because of stronger edge heat- 
ing at higher collision frequencies which prevent more en- 
ergy transported from the antenna into the core plasma. 
The amplitude profile at f e fF = 9.5i/ e j, which agrees with 
experimental data, shows consistent feature with a pre- 
vious study that the RF energy is almost all absorbed 
in the near region of the antenna rather than in the far 
region.— We further studied the effect of the direction 
of static magnetic field on wave propagations, and found 
the antiparallel feature that waves propagate in the op- 
posite direction of magnetic field. This dependence of the 
direction of helicon wave propagations to that of static 
magnetic field in a non-uniform field configuration is con- 
sistent with previous observations made in uniform field 
configurations l 11 - l 22 2Ml 

Physics questions raised by this work include: fur- 



ther explanation of exactly how axially non-uniform field 
might affect the radially localized helicon mode^ inclu- 
sion of different m numbers in the glass layer and any 
subsequent coupling to the plasma at the plasma-glass in- 
terface, and identification of independent first and second 
radial modes that superpose to yield a local minimum in 
wave field amplitude at z = 0.27 m. Experimental mea- 
surements that might corroborate the wave field genera- 
tion mechanism and associated physics include: the mea- 
surement of axial profile of density, measurements of n e , 
B, and E with a reversed field. 
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